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ABSTRACT 

Over two thousand Galactic microlensing events have been discovered so far. 
All of them can be explained by events caused by single or multiple lenses (in- 
cluding binaries and planetary companions). However, when a microlensing event 
occurs in highly dense star fields such as in the Galactic bulge or in a globular 
cluster, it is necessarily affected by the shear from the global distribution of 
mass near the lens star. We investigate the distortions due to this shear in the 
microlensing light curves and in the astrometric microlensing centroid shift tra- 
jectories. As expected, the light curve deviation increases as the shear increases 
and the impact parameter decreases. Although the light curve in the presence 
of a small shear is similar to the simple Paczyhski curve with a slightly smaller 
impact parameter, the detailed difference between the light curve with and with- 
out shear reflects the direction and the magnitude of the shear. The centroid 
shift trajectory also deviates from a simple ellipse in the presence of shear. The 
distortion of the centroid shift trajectory increases as the impact parameter de- 
creases, and the shape of the trajectory becomes complicated when the impact 
parameter becomes small enough. The magnitude of the maximum distortion 
depends on the magnitude and the direction of the shear. For a source trajec- 
tory in a given direction, the time of the maximum distortion depends mostly on 
the impact parameter and hardly on the shear. It is possible to determine the 
magnitude of the shear and its direction if both the time and the magnitude of 
the maximum astrometric distortion are measured. The magnitude of the shear 
produced by the Galactic bulge or a globular cluster falls in the range 10~^-10~^ 
in normalized units. Although the actual determination of the shear from the 
Galactic sub-structures will not be easy due to complications such as binary com- 
panion, future large scale microlensing experiments may enable us to determine 
the shear in some high amplification events, leading eventually to mapping the 
Galactic mass distribution. 
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Subject headings: gravitational lensing — Galaxy: bulge — globular clusters: 
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1. Introduction 

Since Paczynski (1986; see also Gott 1981) first proposed gravitational microlensing as a 
tool to detect massive astronomical compact objects (MACHOs) in the Galactic halo, three 
groups (OGLE, Udalski et al. 1992; EROS, Aubourg et al. 1993; MACHO, Alcock et al. 
1993) independently discovered the first microlensing events, and subsequent observations 
detected by now more than two thousand microlensing events. 

Gravitational microlensing has been applied to various fields of astronomy, such as the 
studies of Galactic structure and stellar populations and the search for extra-solar planets. 
With increasing potential of microlensing as a versatile astrophysical tool, more advanced 
microlensing experiments such as highly precise follow-up observations and pixel lensing 
observations are currently being carried out, and next generation microlensing experiments 
such as astrometric microlensing observations by using the Space Interferometry Mission 
(SIM) and Keck and VLT interferometers have been proposed. 

Microlensing experiments are conducted toward very dense star fields such as the Galac- 
tic bulge and the Magellanic clouds. When a microlensing event occurs in these crowded 
fields, it is ne cessarily affected by the shear caused by the global distribution of mass around 



the lens star. IChang fc Refsdall (119791 . 1 19841 ) have discussed the effect of a star on the macro- 
lensed image produced by the galaxy as a whole, often referred to as the 'Chang & Refsdal 
lens'. Their lens model is characterized by 'convergence' and 'shear': the former depends on 
the mass density within the beam and determines the magnification of the image, while the 
latter depends on the mass distribution ou tside of the beam and determines the distortion 



of the image (see e.g. ISchneider et al.lll992l ). They pointed out that the configurations and 
the observational characteristics of macro-lensed images can be significantly affected by a 
single star. However, their work considered quasar lensing under shear, and focused only 
on the image configurations. In this paper, we investigate how the shear affects the stellar 
microlensing light curves and the astrometric centroid shift trajectories. 

The paper is organized as follows. We first introduce a single point-mass microlensing 
in §2 to compare with the point-mass lensing with shear in §3.1. In §3.2 and §3.3, we discuss 
the distortions in the micro-lensed light curves and centroid shift trajectories in the presence 
of the shear. In §4, we discuss possible applications to Galactic microlensing experiments. 
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In §5, we summarize the results and conclude. 



2. Single Point Mass Microlensing 



Microlensing occurs when a lensing mass passes very close to the line of sight to a 
background source star. When the lensing mass is a single point mass, the source star splits 
into two images. The separation between the two images is of the order of the Einstein 
radius, which is related to the physical parameters of the lens system by 
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where m is the mass of the lens, Di and Dg the distances from the observer to the lens and 
source, respectively, and Dig the distance from the lens to the source. The angular positions 
of the images with respect to the lens are 
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is the projected lens-source separation vector in units of the Be, uq the impact parameter, to 
the time of maximum amplificatio n. The Einst ein ring radius crossing time or the Einstein 
time scale, tE, is given by (see e.g. lGouldll2000l ). 
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Tlrel and firei are the relative source-lens parallax and proper motion, respectively. The 
magnification of each individual image is given by 
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where Ao,+ and Aq_ are the magnification factors of the major and minor images, re- 
spectively. The position of the center of light (centroid) corresponds to the magnification- 
weighted mean of image positions, i.e. 
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with Aq = v4o,+ + Aq,- is the total magnification. The centroid shift S6q is defined as the 
difference between the image centroid 6q and the unlensed source position ^^^o; and is related 
to the lensing parameters by 



56>n = 6>n - 6> 
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(7) 



The position of the centroid shift caus ed by a single point- mass lensing follows an ellipse 
(jWalkerl 1 19951 : iJeong. Han fc ParkI Il999l ) , which is represented by 



^ + 1 = 1. (8) 

where the a; and y represent the centroid shift parallel and normal to the lens-source transverse 
motion, respectively. The semi-major axis a and semi-minor axis b depend on the impact 
parameter uq as 
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3. Microlensing under Shear 
3.1. Lens Equation 



When a source is len sed by a point mass m plus p l anar mass distribu tion, the lens 
equation becomes (see e.g. ISchneider et al.lll992l : lAnll2005l ; lAn &: Evand l2006l ) 



r 



a(r\ 



(10) 



where the two-dimensional vectors r and s are the positions of the images and the unlensed 
source, respectively. The vector r in the lens (image) plane is normalized by te = OeDi and 
the vector s in the source plane by s^; = OeDs- The scaled deflection angle a(r) due to 
the additional mass distribution in the lens plane is given by the gradient of the deflection 
potential ip: 

a{r) = Vipir) , (11) 



where 



d r a(r') In \r — r' 



(12) 
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Here (T{r) represents the surface mass density S(r) normalized by the critical surface mass 
density Scr, 

. X S(r) c'^Ds , , 
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where Dp is the distance from the observer to the planar mass distribution, which is in 
general different from Di, the distance to the point mass m, and similarly for D^/g. 

When the mass distribution consists of a stellar mass and a much larger-scale extended 
mass distribution, the lens equation is approximated by a point mass plus quadrupole lens 
model (Chang & Refsdal 1984; see also Kovner 1987), 



r-^-(^ + ^ ° ]r. (14) 



The quadrupole term is specified by the convergence k and the shear 7. These quantities 
are the two-dimensional second derivatives of ip{r): 

V'll + -022 , ,T .X 

« = 2 ' 7 = V7i + 72 , (15) 

where 

V'll-'022 , , 

7l = ^ , 72 = -012 = -021 , (16) 

and ipij is the partial derivative of il'{r) with respect to r^, 
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If we define the rescaled coordinates S and R as 

S ^ ^ , R = V|l-/^ + 7k , (18) 

V|l -K + 7I 



equation (14) becomes 



^ , A \ „ i^ 



where 

£ = sign(l-«: + 7), A = |^, (20) 

1+7 

and the reduced shear 7 is 

For convenience, 7 will be referred simply as the shear henceforth. 

In order to solve the lens equation, we introduce the polar coordinates {R, ip) in the lens 
plane: = Rcos(f and Ry = Rsinip. Then equation (19) becomes 

= R{eA - R-^) cos ip, Sy = R{£l - R""^) sin p) , (22) 
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which yields a fourth-order equation for (jSchneider et al.lll992l ). 



K^R^ - [£2A(A+ 1^ 

+ [A2 + 4A + 1 + e2{Sl + ^Sl)]R'' 
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We solve equation (23) by Laguerre's method. 
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3.2. Light Curve 

Equation (23) has either zero, two, or four real roots, each of which corresponds to the 
position of the individual image. The Jacobian matrix of equation (19) is 
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whose determinant is 
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The magnification of each individual image is 
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The total magnification is the sum of the magnifications of individual images, = A^^-. 
In order to quantify how much the light curve in the presence of shear deviates from that in 
the absence of shear, we define the excess magnification as 



5A = A^- Ao 



(27) 



where A^ and Aq represent the magnifications with and without shear, respectively. In Figure 
1, we present the contour maps of magnification A^ (left panels) and excess magnification 
6A (right panels) as a function of source position {sx,Sy) for 7 = 10~^, 10""^, and 10~®, 
respectively. The caustics appear as the central diamonds in the left panel of Figure 1, 
whose full width on the x-axis is 47(1 — k + 'j)~^^'^ and that on the y-axis is 47(1 — k — 7)^-*^/^ 
(IHan et al.ll2005l ). When the source is outside the caustic, i.e. uq > 47, the number of images 
is two as in the simple lensing without shear. 
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The series solution of equation (23) can be derived under the assumption k <^ 7 <^ 
\s\ <^ 1. The excess magnification SA is calculated in powers of 7, and the leading term 
yields 



Although other geometric configurations will yield different values, the order of magnitude 
of the deviation will be the same, SA ~ ^/uq. This shows that even very small shear can 
produce a significant deviation in the light curve if the impact parameter is small enough, 
that is in high-magnification events. 

We now investigate the light curves for typical source trajectories. Since the trajectory 
of the source does not coincide with the direction of the shear in general, we choose source 
trajectories with various angles (see Fig. 2). Figure 2 shows the angle ■»? defined as the 
angle between the x-axis and the source trajectory. The dotted ring around the lens (its 
position marked by 'x') is the Einstein ring. In the left panel of Figure 3, we present the 
lensing light curves for the corresponding source trajectories marked in Figure 2 with values 
of the shear 7 = 0.01, 0.005, 0.001, and 7 = 0.0 (no shear). All light curves have the same 
impact parameter, uq = 0.3. The minimum magnification is greater than the minimum 
of Aq by a factor 1/(1 + 7) (1 — 7) for 7 < 1. This is because the brightness of the source is 
increased by the shear alone even in the absence of the point-mass lens. In the right panel 
of Figure 3, we present the magnitude and pattern oiSA for various shear values and source 
trajectories. We find that the excess magnification increases as 7 increases and the excess 
becomes maximum when the source trajectory is parallel to the shear direction (1? = 0°). 
Even if the light curves have the same shear and the impact parameter, the positions and 
heights of the peak deviations vary with i). We find that the value of maximum SA decreases 
as source trajectory becomes perpendicular to the shear direction {'d = 90°). We also find 
that unless the source trajectory is parallel or perpendicular to the shear direction, the 
deviation SA becomes asymmetric in the presence of the shear. 

In real observations of microlensing, however, each light curve will be first fitted by a 
theoretical microlensing light curve such as the Paczynski curve. So we check how the light 
variation induced by the presence of shear deviates from the single mass lensing curve. In 
Paczynski curve, the peak amplification Ap is a simple function of the impact parameter uq, 



We compare the light curve in the presence of shear with the Paczynski curve that has 
the same peak magnification Ap at t — to against the background magnification due to the 





1/2 



(28) 



shear alone, A. 



. The impact parameter Uop (eq. [28]) from the fitted Paczynski curve 
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is always slightly smaller than the true impact parameter uq because shear increases the 
maximum magnification. The detailed shapes of the deviations from the Paczyhski curve 
for various 7 and are shown in Figure 4. As expected from the Figure 3, the light curve 
becomes asymmetric in the presence of the shear when the source trajectory is not parallel 
or perpendicular to the direction of the shear. By fitting the deviation from Paczynski curve 
with an appropriate lens model with shear, it is possible, in principle, to determine the 
magnitude and the direction of the shear if the photometric accuracy is good enough to 
measure the deviation. In real observations, the light curve as a whole will be fitted by the 
Paczyihski curve, and the shape and the maximum value of the deviation can be somewhat 
smaller than that from the simple fitting of the maximum amplification. 



3.3. Centroid Shift 

We define the centroid shift by 

Sd^ = e^-es, (29) 
where ^ 

^■T ^ 'a^ ~ ' ^ ^^^^ ' ^"^^^ 

s-y is the position of the image that would result from the shear alone in the absence of 
the point mass, i.e. = s + cx.{s^). When microlensing occurs in the presence of shear, 
astrometric observations will measure Sd^. From equations (18) and (22), we see that 

*7,a; = ) ^7,2/ = ; ) (31) 

1 — K — 7 1 — K + J 

which shows the whole field is sheared even before stellar microlensing occurs. This affects 
the observed proper motion of the source as well as the estimate of the lensing parameters. 
If the source proper motion increases by the factor A~^, the Einstein ring radius crossing 
time would change to 

t'E ^Me. (32) 

The trajectories of Sdj in units of 9e are shown in Figure 5 for given and 7's. The 
left panel shows the change of centroid shift trajectories depending on ■»? and 7. The centroid 
shift trajectories {dotted, dot-dashed, and dashed curves) deviate from a simple astrometric 
ellipse {sohd curve), and both the shape and the magnitude of the distortion vary with 
and 7. The right panel of Figure 5 shows the variation of the centroid shift trajectory 
depending on the impact parameter uq. The trajectories of centroid shift for different Uq in 
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the presence of the shear are different from each other. The shape becomes more complex 
and the distortion increases as Uq decreases. The distortion becomes very large for small 
values of Uq even when the shear is small. 

To quantify the deviation of the centroid shift trajectory, we calculate the excess centroid 
shift defined as 

Ad = 56^ - SOo , (33) 

where 66^ and 66q represent the centroid shifts with and without shear, respectively. The 
deviation A6 is really the deviation from the centroid shift ellipse which is expected in a 
point mass lensing without shear. The trajectories of A6 are shown in the left panel of 
Figure 6. The arrow in each panel shows the direction of the centroid motion with the 
progress of time. All excess centroi d shifts ha ve one twist. This is similar to the planet- 



induced micro lensing centroid shift (jHaru 120021 ). The magnitude of the excess centroid shift 
defined as A6 = \S0^\ — \66o\ is shown in the right panel of Figure 6. We find that major 
deviation occurs when —2 < (t — to)/tE < 2, as expected. The sign of A6 can be either 
positive or negative. The detailed shape of A9 as a function of (t — to)/tE depends on {}. 
For intermediate values of i), e.g., = 30° or 60°, AO changes its sign near t = to. 

We further investigate the magnitude of the maximum distortion, A6max, and the time 
of maximum distortion, tmax- Figure 7a shows how AOmax in units of 6'^; varies as a function 
of 7 for different values of {} and uq. For small enough value of 7, \ogA9max increases 
linearly as log 7 increases. For a given 7, A6max becomes maximum when = 0° and 90° 
and minimum when = 45°. Figure 7b shows the dependence of A6max on uq for given 7 
and = 45°. Figure 8a and 8b show tmax in units of t^; as a function of 7 for different i) and 
Uq, which Figure 8c shows the dependence of tmax on uq. For each and uq, tmax is nearly 
constant, independent of 7. Since Uq can be determined from the centroid shift trajectory, 
7 can be determined from the AO^ax (Fig. 7a) and ^ from tmax (Fig. 8a). Therefore, it 
is possible to determine the shear and its direction if we determine by astrometry both the 
time tmax and the magnitude of the maximum astrometric distortion A9max- Needless to say, 
one can always fit the full lensing (with shear) model to each individual case, and determine 
the shear and its direction. 

As in §3.2, we can also calculate the series expressions for in powers of 7 under the 
same assumption: 



A6' ~ 7 

where 
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(34) 



+ Ug, Sy = —T sin + Uq COS (35) 



and 
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t_k 

tE 



(36) 



The time of maximum distortion, tmax, is given by a root of the equation 

+ {'iul + 1)T^ + {2ul + Auo) tan 2i? + {3u^ - 5ul)T + {2ul - 2ul) tan 2i? = 0. (37) 

This equation has three real roots, and the one with the smallest absolute value corresponds 
to tmax- The angle i? is now 



— - arctan 
2 



{2ul + Auo)tl,, + {2u^,-2ul) 



(38) 



Hence, we can also approximately determine 7 and its direction from equations (34) and 
(38) from tmax and A6max- 



4. Application to Galactic Microlensing 

Now we discuss Galactic microlensing affected by the shear. Consider a microlensing 
system that consists of a single lensing star under a shear. The shear field can be produced 
by any Galactic sub-structures such as globular clusters and the Galactic bulge. Here, we 
only consider the Galactic bulge and globular clusters as typical examples of the shear, and 
model their mass distribution as a point mass or the Plummer's model. The Galactic bulge 
is significantly extended along the line of sight and, therefore, the mass distribution has to 
be weighted by the factor [DiiDiig)/ Dg and projected along the line of sight; the bulge mass 
distribution located near the source plane contributes little while those near the half of the 
distance to the source contributes most. But in this work we model the bulge as a planar 
mass distribution at the same distance as that of the lensing star for simplicity. 

1. Point mass: When the lensed images are located far from the source of the shear, 
the shear field may be approximated by that produced by a point mass. When the mass 
distribution consists of a lens with mass m located at the origin of the coordinate and an 
additional mass with mass M located at ri (in units of r^), the (additional) deflection 
potential '4>{r) in equation (12) becomes 

M 

ip(r) = — Inlr — ri|. (39) 
m 

Shear 7 and convergence k are calculated from equation (15): 
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Convergence disappears (k = 0) because the beam is empty, and the reduced shear is equal 
to the shear, 7 = 7. 

2. Plummer's model: As an example of extended mass distribution, we choose the Plum- 



mer's model that approximates the surface mass distribution of a globular cluster (jPlummer 



19151 : iBinney fc Tremaind 119871 ) . In the Plummer's model, the surface mass density is ex- 



pressed as 

..Ml 

S(r) = — ^- -— — -TTTTT , (41) 

^ ' vrrg [1 + (r/ro)2]2 ' ^ ' 

where M is the total mass and tq is the core length in units of rg. The radius rh containing 
half the total mass is equal to 1.3048ro. Then, the deflection potential V'l^) due to the 
Plummer's mass distribution centered at becomes 

ip{r) = — / dr — ^ -m\r — r\ 

= ^\n[{r-r^f + TlY/K (42) 
Again, 7 and /t are calculated from equation (15): 



M 

m (ri^ + Tq)^ 



7 = -TTT^ (43) 



and 

- = -T-fhrl ■ (44) 

Unlike the point mass case, the convergence k is not zero because there is mass within the 
beam. 

Figure 9a shows the values of 7 and /t as functions of the distance from the center of 
the mass distribution when the shear is produced by the Galactic bulge with a total mass 
of M = 1.3 X 10^° M0 and located at 8.5 kpc from the Earth. We assume that the lens 
star of IMq is also located at 8.5 kpc and the source star at 9.5 kpc. Solid curve represents 
7 when the bulge is modeled as a point mass while dotted curve as the Plummer's model 
with ro = 500 pc. Dashed curve shows k for the Plummer's model. Figure 9b is the same 
as Figure 9a for a 10^ M© globular cluster at 4 kpc with tq = 2 pc, 1 lens star at the 
same 4 kpc, and a source star at 8.5 kpc. Hence, we expect the shear 7 produced by typical 
globular clusters or the Galactic bulge to be in the range of 10~^ ~ 10~^. 

Shown in Figure 10a is the maximum distortion of the centroid shift in arcseconds when 
a microlensing event occurs in the Galactic bulge. If future astrometric observation can 
achieve the positional accuracy down to ~ 1 micro-arcsec, then 7 ~ 10~^'^ shear field can 
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be detected in very high magnification events with uo < 0.002. When a microlensing event 
occurs near a globular cluster (Fig. 10b), 7 ~ lO"'''^ shear field can be detected for events 
with Uo < 0.05. On the other hand, in order to detect the deviation of a typical light curve 
for 7 — 10~^-^ and Uq — 0.01, photometric accuracy should be better than Am < 10~^-^. 

Although we may be able to measure the shear in microlensing events, the measurement 
alone docs not tell us about the source of the shear. The shear may be from the Galactic 
sub-structures, but it can also be from many other objects. The most obvious source is the 
binary companion. A binary companion at a typical distance of ~ 35 AU (Duquennoy & 
Mayor 1991) will produce a shear close to 10^^, two to three orders of magnitude larger than 
the shear by the Galactic bulge. The distribution of the binary period and the mass ratio 
(Duquennoy & Mayor 1991) implies 85% of all binary companions will produce 7 > 10"^. So 
the shear due to binary companion will dominate or be comparable to the shear expected from 
Galactic sub-structures. The shear can be also produced by planetary companions. Planets 
detected in the current microlensing experiments produce shear in the range 7 ~ 10~^ — 10~^. 
A typical Earth-mass planet at a distance of ~ 2Re {Re — 4 AU for 1 Mq) from the lens 
produces 7 ~ 10^^, while Pluto-mass planet at a distance ~ lOi?^; produces 7 ~ lO^^''. 
Projected companion stars, unassociated but located near the lensing star in the projected 
sky, will also affect lensing similarly. Hence, in practice it will be difficult to identify the 
shear by Galactic sub-structures against the shear by companion stars. Still, there can be 
some microlensing events that are not affected by companion star or in which the shear by 
Galactic sub-structure may be measured from the statistical analysis of many events as in 
cosmological weak lensing systems. 

There may also be cases for which we may set an upper limit on the value of the shear. 
A very small upper limit on shear suggests the non-existence of a binary companion or a 
Galactic structure. For example, if we measure 7 < 10~^-^, we can expect that the lens does 
not have a companion with more than 95% confidence if we assume the distribution of the 
binary companion by Duquennoy & Mayor (1991). 

Complications in usual microlensing events, for example, blending and binary source, 
will also affect microlensing under shear. For example, blending decreases the deviation in 
the image centroid shift and the light curve (Fig. lib as compared to Fig. 11a) and the 
binary source completely messes up the centroid shift trajectory (Fig. 11c). However, since 
the shape of the trajectory mainly depends on the size of the shear, fitting the full trajectory 
may sort our these complications. 

Although the typical shear expected from the Galactic bulge or globular clusters is 
too small or the clear case that enables the determination of shear is too infrequent to be 
comfortably detected by current or near-future microlensing experiments, next generation 
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microlensing experiments may enable us to measure the magnitude and the direction of 
even smaller shear fields among numerous microlensing events. Then, from the theoretical 
framewo rk of weak lerising, the shear map can be inverted to reproduc e the mass distribution 



(see e.g.. lMellierlll999l : iBartelmann fc Schneider! l200ll : iRefregierl |2003| ) . making it possible to 
map or at least constrain the Galactic mass distribution. Since all our discussions are based 
on normalized units, the reproduced mass distribution will be in units of Ecr per 6^. If Dg 
and Di are additionally determined, the mass distribution can be determined in physical 
units. 



5. Summary 

We investigated microlensing under a shear, which might be produced by Galactic sub- 
structures such as globular clusters or the Galactic bulge. We analyzed its effect on the 
microlensing light curves and astrometric centroid shift trajectories. We found the foUowings: 

1. The light curve deviation from the Paczyhski curve increases as the shear increases 
and the impact parameter decreases. The positions and heights of the maximum deviation 
vary depending on the direction of the source trajectory: the light curve becomes asymmetric 
if the source trajectory is not parallel or perpendicular to the shear direction. 

2. The centroid shift trajectory in the presence of shear deviates from a simple ellipse, 
especially when the source is within the Einstein ring. The magnitude of the maximum 
distortion depends on shear and its direction, and becomes largest when the trajectory is 
parallel or perpendicular to the shear direction. The time of maximum distortion is nearly 
independent of the impact parameter. The distortion of the centroid shift trajectory increases 
as the impact parameter decreases and the shape becomes very complex when the impact 
parameter is very small. 

3. If we measure the distortion of the astrometric centroid shift trajectory and the time 
of the maximum distortion, we can determine the shear and its direction. 

4. The magnitude of the shear produced by the Galactic bulge or globular clusters near 
the Galactic center is of the order of 10~^ to 10~^ in normalized units. This shear, in principle, 
could be detected by future microlensing experiments, especially in high magnification events. 

Successful measurement of the shear in various directions in the Galaxy with next gen- 
eration microlensing experiments may eventually lead to the mapping of the Galactic mass 
distribution. 
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Fig. 1. — Magnification map. Grayscale maps of tiie magnification (left panel) and the 
excess magnification 5A (right panel) as a function of source position {s^, Sy) for reduced 
shear 7 = 10~^ (top), 10~^ (middle), and 10~^ (bottom). Contours in the left panels represent 
Aj — 40 (long dashed curve) and Aj — 80 (dotted curve). Gray scale in the right panels 
represents positive (bright) and negative (dark) deviation regions and dot-dashed curve shows 
SA — regions. 
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Fig. 2. — Geometry of lensing. The straight hnes with an arrow represent source trajectories 
with different ■&. We define ■& as the angle between the a;- axis and the source trajectory. All 
source trajectories have the same impact parameter uq — 0.3. The dotted ring around the 
lens, marked by 'x', shows the Einstein ring. 
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Fig. 3. — Light curves and the excess magnifications. The left panels show the light curves 
for the source trajectories {uq = 0.3) in Figure 1 with the shear value of 7 = 0.01, 0.005, 0.001. 
All curves show the magnification near the peak, and the solid curves are the magnification 
without shear (7 = 0.0). The right panels show the difference in magnification between the 
light curve in the presence of shear and that in the absence of shear as a function of time 
for three different source trajectories. 
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Fig. 4. — Deviation in the light curves. Difference between the hght curve with shear and the 
corresponding Paczyhslci curve as a function of time for different source trajectories. Lines 
are for the same 7's in Figure 3. 
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Fig. 5. — Centroid shift trajectories. The left panels show the change of centroid shift 
trajectories in units of for different angle and shear 7, while the right panels show the 
variation of centroid shift trajectories for different impact parameter uq. The solid curves 
represent the simple astrometric ellipses. 
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Fig. 6. — The trajectories and magnitudes of excess centroid shifts. The left panels show 
the trajectories of excess centroid shifts of the events for 7 = 0.01 cases in the left panels 
of Figure 5. The arrow mark shows the direction of event progress. The right panels show 
the magnitudes of excess centroid shifts in the left panels as a function of time with three 
different values of shear, 7 = 0.01, 0.005, 0.001. 
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Fig. 7. — The maximum distortion of excess centroid shift, (a) The maximum distortion 
of excess centroid shift ^9 max for different values of Uq and ■»? as a function of shear 7. We 
show I^Omax for ■& — Q° and 45° only because other source trajectories appear between the 
two cases, (b) The maximum distortion i^^Omax as a function of uq for different values of 7 
and 1? = 45°. 
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Fig. 8. — The time of maximum astrometric distortion. The time of maximum distortion 
of excess centroid shift, tmax, as a function of 7 and uq: (a) for fixed uq — 0.3, (b) for fixed 
1? = 30° and (c) for fixed 7 = 10-^. 
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Fig. 9. — Shear and convergence from the Galactic sub-structures, (a) Shear 7 and conver- 
gence K, as functions of the distance from the center of the mass distribution for 1.3 x 10^° 
M0 Galactic bulge at 8.5 kpc for 1 Mq lens at 8.5 kpc and a source star at 9.5 kpc. Sohd 
curve shows the value of 7 when the bulge is modeled as a point mass while dotted curve as 
a Plummer's model with tq = 500 pc. Dashed curve shows k, for Plummer's model, (b) The 
same for a 10^ Mq globular cluster at 4 kpc with ro = 2 pc and a source star at 8.5 kpc. 
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Fig. 10. — The maximum astrometric distortion expected. The maximum distortion AOmax 
in arcseconds in excess centroid shift as a function of the shear 7 due to (a) the Galactic 
bulge {Dg = 9.5 kpc, Di = 8.5 kpc) and (b) a globular cluster {Dg = 8.5 kpc, Di = 4.0 kpc), 
both with a 1 Mq lens. Solid lines are for -& — 0°, dashed ones for ■»? = 30°, and long dashed 
ones for 1? = 45°. 
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Fig. 11. — Complications in the ccntroid shift trajectory and the hght curve from blending. 
The centroid shift trajectories (left), centroid shift deviations (middle), and light curves 
(right) are presented in case of (a) no blending (top), (b) bright lens (middle), and (c) 
binary source (bottom). The light fraction of the blending source is 0.5 for the lensing 
parameters Uq — 0.1, ■»? = 30°, and 7 = 0.01. 



